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This is the first of two papers on computing the self-force in a radiation gauge for a particle of 
mass m moving in circular, equatorial orbit about a Kerr black hole. In the EMRI (extreme-mass- 
ratio inspiral) framework, with mode-sum renormalization, we compute the renormalized value of 
the quantity H := ^h a pu a u^ , gauge-invariant under gauge transformations generated by a helically 
symmetric gauge vector; here h a p is the metric perturbation, u a the particle's 4-velocity. We find 
the related order m correction to the particle's angular velocity at fixed renormalized redshift (and 
to its redshift at fixed angular velocity), each of which can be written in terms of H. The radiative 
part of the metric perturbation is constructed from a Hertz potential that is extracted from the 
O f Weyl scalar by an algebraic inversion We then write the spin- weighted spheroidal harmonics 

D ' as a sum over spin- weighted spherical harmonics s Yi m and use mode-sum renormalization to find 

^0 ' the renormalization coefficients by matching a series inL = ^+ l/2to the large- L behavior of the 

expression for H . The non-radiative parts of the perturbed metric associated with changes in mass 
and angular momentum are calculated in the Kerr gauge. 
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I. INTRODUCTION 



We report here a first computation of the order m renormalized corrections Af2 to the angular velocity of a mass 
[ m in circular orbit about a Kerr black hole and of the related change AU in the scalar U — u a V a t. This is the 
0^ ■ constant of proportionality between the particle's velocity and the helical Killing vector k a tangent to its trajectory: 
' u a = Uk a . The corresponding computation for a Schwarzschild background was first done by Detweiler [2|, who 
. pointed out the invariance of these quantities under gauge transformations generated by helically symmetric gauge 
f*"*- " vectors. The invariance allows straightforward comparison with computations in other gauges, and comparisons were 
done for circular orbits in Schwarzschild for a Lorenz gauge [|| and for our modified radiation gauge [1] . Previous 
conservative self-force computations in Kerr spacetime have been restricted to a scalar particle in eccentric, equatorial 
orbit 

The radiation gauges are associated with the Weyl scalars ipo and ■04 that satisfy the Teukolsky equation Q : They 
are gauges in which one can construct a perturbed metric in the vacuum spacetime outside the particle from ipo 
or ?/>4, and for circular orbits this can be done by an algebraic inversion. The resulting metric, however, lacks the 
non-radiative contributions from the changes SM and SJ in mass and angular momentum, and from the change in 
the center of mass. We write the metric perturbations corresponding to SM and 5 J in the Kerr gauge, the result of 
changing a and M in the Kerr metric written in Boycr-Lindquist coordinates. It is in this modified radiation gauge 
that we compute Af2 and AU. 

We believe that the numerical results given here are accurate to better than one part in 10 8 , consistent with 
our earlier Schwarzschild computation. Numerical values we reported in an earlier eprint version of this paper had 
significant inaccuracy, due to an oversight in our computation of the singular field, found after two comparisons. We 
performed a detailed comparison of our results for H rcn with corresponding post-Newtonian values found by A. Le 
Tiec; and we compared a few values of iF cn with those of S. Dolan, who has completed another EMRI calculation 
of H rcn . Dolan's calculation uses a Lorenz gauge and an effective-source renormalization method (see, e.g. [cIlTil]). 
Where our initial error led to a numerical error of about 5% (with about twice that error for the the difference between 
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H for Schwarzschild and for Kerr), Dolan's present accuracy is about 1%. 1 

The plan of the paper is as follows. In Sec. Ill Al after a brief review of the Teukolsky equation and associated 
formalism, we obtain the analytic form of the Weyl scalar ^>5 et , expressed as a sum over spin- weighted spheroidal 
harmonics. In III Bl we review the construction of the metric perturbation in a radiation gauge, using a Hertz potential 
constructed algebraically from tj)Q. We describe the mode sum renormalization, pointing out a difference, not present in 
our Schwarzschild computation, between the singular field at subleading order in a Lorenz gauge and in our radiation 
gauge. In Sec. Ill CI we first review two different series expansions for spin- weighted spheroidal harmonics, one as 
a sum of Jacobi polynomials, the other as a sum over spin-weighted spherical harmonics. We check the accuracy of 
the angular harmonics by computing them using these two different expansions. We then obtain explicit expressions 
for the values of the harmonics {Ye m and & t 8 = k/2 (the plane of the particle's orbit). In Sec. HID I we 

review the numerical construction of solutions to the radial Teukolsky equation, by direct numerical integration and 
by numerical integration of the Sasaki-Nakamura equation. Again we check the accuracy of our computation of the 
radial harmonics and the angular harmonics of ipo ct obtained in these two independent ways. In Sec. Hlll we compute 
the tetrad components of the perturbed metric from the Hertz potential. In Sec. II VI we find expressions for the change 
AU in the renormalized redshift factor at fixed fl and the corresponding change in the angular velocity at fixed U ; 

each is simply related to the quantity H lcn = -h Y °pU a u^ . In Sec. |V]we present the computation of the retarded 

component of the lower multipoles in Boyer-Lindquist coordinates. In Sec. IVII we present the numerical results for 
H ren , AU and Afi. Finally, in Sec. VII, we briefly discuss our results and future work. 



II. REVIEW AND COMPUTATION OF 4> { 
A. Formalism 



We work in Boyer-Lindquist coordinates where the Kerr metric is given by 

/ 2Mr\ , , AMarsm 2 9 , ,, £,„ ~ / 9 , 2Mra 2 sin 2 6>\ 2 „ , ,o 
ds 2 = I 1 — J dt 2 + dtdcj) - —dr 2 - Y>d0 2 - [r 2 + a 2 + J sin 2 6d(j> 2 , (1) 

with A = r 2 — 2Mr + a 2 and £ = r 2 + a 2 cos 2 8. The Kinnersley tetrad vectors have components 

r = i(r 2 + a 2 ,A,0,a), 
n Q = ^(r 2 +a 2 ,-A,0,a), 

m a = ^§ (ia sin 0,0,1, -M , (2) 

and we denote by D,A and 8 the derivative operators along the tetrad vectors l a ,n a and m a , respectively. The 
non-vanishing spin coefficients associated with this tetrad are 

— 1 —gcot9 iag 2 sin —ias'mO Ag (r — M) - 

r — iacosO 2\/2 y2 v2E 2S 2S 

(3) 

Consider a particle of mass m orbiting a Kerr black hole in circular, equatorial orbit with radial coordinate r = tq. 
As noted in the introduction, the particle's velocity is tangent to a helical Killing vector, 

u a =u t k a , (4) 

where 

k a = t a + Q,cj) a , (5) 



Dolan's results are to be reported in a paper in preparation coauthored by L. Barack and B. Wardell. 
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with 4> a and t a the rotational and asymptotically timclike Killing vectors of the Kerr geometry and Q the particle's 
angular velocity measured by an observer at infinity. For a circular geodesic, the values of u and are 



r 3/2 ± M X/2 C 



n 



±M 1 / 2 



(6) 



rl /2 ± oM V2 ' 



where the upper (lower) sign correspond to direct (retrograde) orbits and the corresponding stress-energy tensor is 
given by 



m 



L a u 5{r - r )6(8 - Tr/2)5((f> - fit). 



(7) 



The perturbed spin-2 Weyl scalar, ipo — —C a i3 1 sl a mPPm s , satisfies the Teukolsky equation 



(r 2 +a 2 ) 



2 i2 



2 • 2 

a" sin i 



9Vo 4A/ar 9 2 ^o 



a(r — M) icos( 



sin 2 9 



dt 2 
d(f> 



A dtd(j) 



1 



sin 2 e 



d 2 ip _ _t__S 
deb 2 ~ A 2 ~dr 



+ 2(2cot 2 6»-l)7/» -4 



M(r 2 - a 2 ) 



A 



r — zacosl 



A 3^o 
dr 

dtpo 
dt 



_J_d_ 

^9 89 

8ttET, 



sin 



dtpo 



(8) 



whose source T is given by 



T := - («5 + w - a - 3/3 - At) {5 + w - 2a - 2/3) Tn 

+ [(5 + w - a - 3/3 - At) (D - 2g) + (D - Ag - g) (S + 2m - 2/3)] T 13 
-{D-Aq-q)(D- q) T 33 

_ T (0) +T (1) +T (2) 



(9) 



where the T^ 1 ' label the three terms in the expression for T . 

Because the Teukolsky equation is separable in the frequency domain, we can write tpo as a sum of angular harmonics, 



-00 = X] R i"i( r ) 2Si m {0)e 



im(<p — Q,t) 



(10) 



where we used the fact that, for circular orbits, harmonics e lm ^ have frequency u = mil. For r =^ ro, the source 
vanishes, and R( m satisfies the radial equation 



AR'; m + 6(r-M)R' h 



K 2 - AiK(r - M) 



+ 8imilr — 7 + 2mj + 6 — 2Eg m I Rt 



0. 



The function iSg m satisfies the angular equation 



1 d 
sm~6>d6> 



dO 



9 9 Am cos 9 + A + m 2 
Y cos z 9 - 4 7 cos 9 h 2 E tm } 2 S em = 0, 



sm 



where A' = mfl(r 2 + a 2 ) — am and 7 = amfl. 

The tetrad components of the stress-energy tensor that enter the expression for T in Eq. (|5J| are 

Tu = OT(1 ~^ - ^o)<5(^ - tt/2)5(0 - fit), 



r 3 



im(l - afi) (a - fi(r^ + a 2 )) u* 



V2rg 

(a-ft(r 2 +a 2 ))V 



S(r-r )5(e-Tr/2)6(cf)-nt), 



S(r-r )S(e-'K/2)5(<l>-nt). 



(11) 



(12) 



(13) 
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Using Eqs. <j2j) to write the derivative operators and Eqs. ([3]) for the spin coefficients, we obtain explicit forms for the 
source terms TM of Eq. ©. Defining the quantities 



m a — ii r + a 



')] 



qo := q\r=r 



m [a ~ fl (r 2 + a 2 )] 



and writing 



oo 1 



} irn((fi— Q,t) 



(14) 



(15) 



we find 



i °° -2 

T(0) = __L m ^(l-an) 2 u t e lm ^ nt Hs"(0-7r/2) 



47T '■ — ' 7 

m— — oo 



mattfl + sinf9) — m — cot 8 H 4iagsin(9 — iapsin( 

sin tq 



77;,(1 — afl) — i- 



ro 



7n(l — aQ) — 5i 



a 



ro 



S(6-Tv/2)}S(r-r ), 



5'(6-ir/2) 



(16) 



T (1 ^--^m V [a-n(r 2 + a 2 )] u*^*-™* { 26' {r - r )S' {6 - tt/2) 

47T TV) I 



[ — - 2g + iq ) S(r - r )S'(9 - tt/2) + ( 2a777^ - 2m + — - — ) <5'(r - r )<5(6> - tt/2) 
r J \ r r 



+ 



^amJ] — 777 



Aia 



amfl — 777 



2ia 



ro / vn> 



5(r-r )6(6-ir/2)}, (17) 



T (2) = __L m 4 [«-^( ? o+a 2 )] 2 ^ m(0 - at) 



47T 1 — ' Tn 

m— — 00 



-S"(r - ro) 



+ 9o ) + — + - \ S'(r - ro) - ( — + iqo J ( — + »9o ) £(r - ro) 
ro r I V 7 *o / V 7 *o 



5(6 -tt/2). (18) 



To compute V>5 et , it is useful to define a Green's function Gg m (r,r') as a solution to 

~ K 2 - MK(r - M) 



AGL. + 6(r - M)G' tm + 



+ 8imQr — 7" + 2777.7 + 6 — 



G, m = 6(r - r'), (19) 



namely 



G im (r,r') =A em A' R H (r < )i? 00 (r > ), 
where A/*, 



A^RhR'^-RooR'h)' 



(20) 
(21) 



Here and i?oo are two linearly independent solutions to the radial equation that are, respectively, ingoing at the 
future horizon and outgoing at future null infinity, and we have suppressed the ^777 indices on Rh and i?oo- With Ga m 
so defined, the full Green's function that satisfies the Teukolsky equation with source 5(r — r')S(cos 9 — cos 8')5(4>— </>') 
is given by 



G(x,x') =Y, A t™A' 2 R H {r < )Roc{r > )2Si m {8)2Si m {6 l )e 



im((j>—(l> ) 



(22) 



tin 
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where we denote by x the spatial point with coordinates r,6,(f> and where A' = r' — 2Mr' + a 2 . Hence, the Weyl 
scalar that satisfies the Teukolsky equation has the form 



ip = -8tt / Y!T{x' ,x )G{x,x')dr'd(cos8')d(f>' 

= -8tt J S'(T (0) +T (1) + T {2) )G{x,x')dr'd{cosO')dcp' 

with £' = r' + a 2 cos 2 6' . The have, for r ^ ro, the explicit forms, 



(23) 



(0) _ 47rm(l - afi^Agu* 

^0 — ~2 



(in 



s , + 2iaS^_ 2So + 2m{an 1)s , + 2iam(an-l)S + ^ 



ro 



(24) 



m 87rim(l — ail) [a — f2(r 2 + a 2 )] u* 



2r 3 



im(tf>—Qt) 



Cm 



2A 2 



roS'n — iaSo + r^So [ amCl — m H 



i?H(ro)i?oo(r)e(r - r ) + i?^(r )i? H (r)e(r - r)] 
+ |2A [S' {5r 2 - 6Mr + a 2 ) - AiaS Q {r ~ M)] - 2A 2 (r S' a - iaS ) 
- S r Al 



im [a - (r 2 + a 2 )] 3 
A r 



Aia\ ( im\a-n(rl+a 2 )] 2 
ami I — m H 1 



y y A 

[a-n(r 2 +a 2 )] 4 



A 



amfi — m — 



ro ) 
2ia 



ro J \ r 

,2 ^ 



+2A S ( amCl - m + - ) (5r 2 - 6Mr + a 2 ) - \R H (r<)i?oo(r>) 

r 2 i 47rmM* (a — £l(r 2 + a 2 )) 2 ^-^ , , , , x . n ,s 

4 = 4^ — Y, A im %S im {e) 2 ft m (7r/2)e™^- n ') x 

r 2 A 2 {i^(n,)i2oo(r)e(r - ro) + iC(r )i? ff (r)e(ro - r)} 



(25) 



4r A (3r 2 - 4Mr + a 2 ) - 2r 2 A 2 



im [a - n (r 2 + a 2 )] 3 
A ro 



{R' H {ro)Roo{r)®{r - r ) + RU^Rh (r)6(r - r)} 



r^A 2 



-r 2 A 2 



im [a — f2 (r 2 + a 2 )] 5 
ro 



Mi 



im [a — (r 2 + a 2 )] 1 
A ro 



-5 2imfir 2zm(r - M)(a - fi(r 2 + a 2 )) 

A 2 

im [a - ft (r 2 + a 2 )] 3 
A r 



§ A 
4r A (3r 2 -4Afr + a 2 ) 



+2(3r 2 - 4Mr + a 2 ) 2 + 4r A (3r - 2M) }R H (r<) J R 00 (r>) 



(26) 



where A = — 27l/r + a 2 , So = 2<%m(7i"/2); Sq and 5 ' are derivatives of with respect to cos0 evaluated at 

e = tt/2. 
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B. Metric perturbation and H len in a radiation gauge 



Before presenting the detailed calculation of H , AU and Af2, it will be helpful briefly to recall parts of the radiation- 
gauge method that we will use. In the CCK procedure [l2rfl6j a perturbed vacuum metric is constructed from a spin 
2 Hertz potential if? that satisfies 

x/j = i [C^ + 12Md t *] , (27) 

where £ 4 = C\LqC-iL-% with C s = 3 S — iasindd t and 3 S = — [dg + icscOd^ — scot 6]. Using Eq (fTUj) and Eq (40) 
of we obtain 

.t, o {-l) m DRl?rn + \2imMVLRl, m im U_ nt) , , , , 

* = 8 2. D 2 + 144Af 2 m 2 tt 2 6 } 2S lm (0), (28) 

where!? 2 = X 2 CH (A c # +2] 2 +8aw( m— aw)Acif(5Aci/+6)+48a 2 u; 2 [2Aci/+3(m— aw) 2 ] and Ac//, the angular eigenvalue 
used by Chandrasekhar [17f , is related to the separation constant iEi m of Eq. (fTTj) by Ac# = + 7 2 ~ 2to7 — 2. 

The perturbed metric is then given (up to parts for which ipo vanishes) by 

h a p = g~ i { n a np(6 — 3a — (3 + 5ir)(5 — 4a + n) + rfi a mp{ A + 5^ — 37 + j)( A + p — 47) 

-n( a mp) [(6 - 3a + (3 + 5n + f)(A + p - 47) + (A + 5p - p - 37 - "y)(6 - 4a + tt)] + c.c. (29) 



Whiting [Til decomposition of the metric in a Lorcnz gauge. This has the form fi°p ho1 = + h™g' LoT , with 

L ° r a smooth solution to the vacuum Einstein equation in a neighborhood of the particle. Because the perturbed 
Weyl scalar ipo is gauge invariant, when computed from this decomposed form of h^t' 1 " 01 , it has the gauge-invariant 
decomposition -0Q Ct — ip^ + ipo en , with ?/>5 en a smooth sourcefree solution to the s = 2 Teukolsky equation. There is 
then a smooth Hertz potential ty ven satisfying the sourcefree s = 2 Teukolsky equation with h™g' given in terms 
of 4 ,ren by Eq. (|2"9"|) . Here /i^™' ORG is the reconstructed perturbed metric up to terms involving changes in the mass, 
angular momentum, and center of mass. (These are metric perturbations for which ipo vanishes, and they correspond 
to the £ = and £ = 1 parts of the perturbation for a Schwarzschild background.) Adding these perturbations in an 
arbitrary gauge that is smoothly related to the Lorenz gauge yields the full perturbed metric in a modified radiation 
gauge. 

We will use this general description at the end of Sec. IIII1 but in rcnormalizing _ff rct = ^h^u a u^ , we exploit 
its invariance under gauge transformations generated by helically symmetric gauge vectors. As in our Schwarzschild 
paper Q , this allows us to use the generic Lorenz-gauge singular behavior of the metric perturbation for the leading 
term in the singular field. A difference between the singular field in our radiation-gauge and in the Lorenz-gauge 
arises at subleading order in the angular harmonic index £ from a gauge vector that is singular at the position of 
the particle, and this is discussed below. The mode-sum renormalization of the metric is described, for example, in 
Sec. IV of a and we briefly review its relevant features. (It is based on the Barack-Ori mode-sum version of the 
MiSaTaQuWa renormalization pjl, [20[ that is reviewed in detail in |2l|.) The components h Y ^' hm of the perturbed 
retarded metric along an ortho normal frame have the singular behavior of a Coulomb field, proportional to p —1 , 
where p is the geodesic distance to the particle trajectory. This implies that iJ rct . Lor has the same Coulomb singular 
behavior. The angular harmonics H 1 g^° 1 Yim of iJ rct > Lor then have finite limits Hg m = lim r ^ ro -ff^' L ° r (£ = 0, r) on 
a sphere through the trajectory of the particle. At the position of the particle, the projection 

ff ret,Lor = £ 0) (30) 

m 

onto the £th subspace has the form 

fl-ret^or = H^ hm + H} en , (31) 

where i?™ n , the value of the renormalized field, falls off faster than any power of £. The singular field H S ' hm ' has 
£- dependence 

#f' Lor = E Q L° + 0(L- 2 ), (32) 
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with Eq independent of L; and the sum over I of the 0(L 2 ) terms vanishes. That is, H s 
Hf = E . The value of H lcn is then given by 



H s + 0(p), where 



— i. r^~; ' » t> m — \ m ' ■■ -» 

=0 



1=0 



S ORG 

In our radiation gauge, we will find that, although H t ' again has the form 



m 



S.ORG 



E L° + O{L- 2 ), 



(33) 



(34) 



with Eq agreeing with its value for a Lorenz gauge, the 0(L 2 ) part gives a contribution of order p° that is odd 
under parity in a hypersurface orthogonal to the trajectory: In particular, the sum of the O(p ) contributions from 



the limits 



and r 



vanishes. 



The loss of gauge invariance at subleading order follows from the way one proves invariance of H ren for a helically 
symmetric gauge vector £ a that is differcntiable at the position of the particle (implying, in particular, iJ ron ^ Lor — 
_ff ron ' ORG ). Under a gauge transformation with u a fixed, H changes by ^£^g a f;u a u 13 = V a £,0U a u 13 . From Eqs. (@| 
and {5J for u a , we have 



(35) 



where we have used the geodesic equation in the last equality. Now u a , as defined by Eq. ((4J, satisfies the geodesic 
equation only on the particle's trajectory; for points a geodesic distance p from the trajectory, it satisfies 



O(p). 



(36) 



The gauge vector relating a Lorenz gauge to our radiation gauge diverges as logp near p = [22|, leading to a 
term of order p° in Eq. (35). We defer to Sec. IIV1 in which the explicit analytic construction of the perturbed metric 
is presented, the discussion of the parity of this term. 

In the actual computation, as described in Sec. (|VI[) . we use spin- weighted spherical harmonics instead of ordinary 
spherical harmonics for parts of the metric with different spin weights, and we then check that the resulting leading 
term in the singular field - the value of Eq - agrees with its value computed analytically within a Lorenz-gaugc 
framework by Linz [23| . 



C. Numerical methods - Angular harmonics 



This section describes two different scries expressions for spin-weighted spheroidal harmonics that we use for two 
different purposes: First, following Fackerell and Crossman 24|, we write each spin- weighted spheroidal harmonic as 
a sum of Jacobi polynomials. The formalism provides an accurate way to evaluate the spheroidal harmonics and the 
angular eigenvalues on which both the angular harmonics and the radial functions Re m depend. Second, to renormalizc 
the metric and the self-force, we write each spin-weighted spheroidal harmonic as sum of spin-weighted spherical 
harmonics. Although the present paper uses only harmonics s S^ m := s Se m u\u= m n, the formalism is developed for 
sStmun with no restriction on uj. 

For 7 = aw, the spin- weighted spheroidal harmonics satisfy the eigenvalue equation 



1 d ( ■ a d a 

( sm9 d~e sSem " 



7 cos 6 — 2s7 cos ( 



2ms cos 9 + s 



sin 



S E(. 



3 SV 



0, 



(37) 



for a Teukolsky-cquation source with time-dependence e lult . The eigenvalue s Ei muj is a continuous function of 7 that 
takes the value £(£+ 1) when 7 = 0. (Note that the spin- weighted spheroidal harmonics and their eigenvalues depend 
on oj only via 7, so the conventional use of the index u instead of 7 is slightly misleading.) 

For fixed s and 7, with 7 real, the eigenfunctions are complete and orthogonal on the sphere, satisfying 



5{cos6 - cos6')8{4> - 4>') = ^ a S tmu {0)e im * s S lm u{0')e- m *\ 

('.711 



(38) 



Jo 



s ^irnuj 



{ey™* s S M (0)e im ' ,t ' sm6d6d4> = 5 U ,8„ 



(39) 
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When 7 = 0, s Simuj{S)e lm ^ becomes the spin-weighted spherical harmonic s Yi m {9, </>), given by 

Y = f [(£-s)\/(£ + s)t /2 B s Y tm , 0<s<£, , . 

™ \ (-l) s l(t + s)\/(£-s)f /2 5- s Ytm, -t<s<0, 

with 

9?7 = — (dg + i esc 6d<f, — s cot 9) 77, 

B?7 = — (dg — i esc + s cot 0) 77 (41) 
where 3 and S arc, respectively, raising and lowering operators for the spin- weight, and rj is a quantity of spin-weight s. 



1. Spin-weighted spheroidal harmonics as a sum over Jacobi Polynomials 



The formalism that expresses the angular harmonic s Se m u in terms of Jacobi polynomials involves several constants 
that depend on parameters s, £, m, and 7 whose values are fixed in this section. To avoid encumbering a large number 
of symbols with the four indices s,m,^,w, in this and the next subsections (Sects. Ill C H and III C 2[) we suppress the 
indices, so that, for example the angular harmonic and its corresponding eigenvalue will be written as 



S S S : 



E := S E£. 



(42) 



To calculate the spin- weighted spheroidal harmonics, their derivatives and eigenvalue, we set x := cos 6* and write 
the homogenous angular equation, Eq. (|37j) . as follows: 



d 2 S 



dS 



(1 - x 2 ) — - 2x— + j 2 x 2 - 2-fsx - 
dx z 



dx 



i 2 + s 2 + 2msx 
1-x 2 



E )S = 0, 



(43) 



The eigenfunctions' dominant behavior at x = ±1 is (1 =F x ym±s\/2^ Following the formalism (and notation) of 
Fackcrcll and Crossman [24]], we introduce a = \m + s\ and /? = \m — s\ for simplicity and introduce new functions U 
and V. as follows: 



S(x) 
S(x) 



ix 



l-X 

2 

1-x 



1 + X 

2 

1 + x 



U(x), 
V(x) 



(44) 
(45) 



The functions U and V satisfy the following differential equation (where the upper sign is used if F = U, the lower 
sign if F = V): 



(1 - x 2 )d 2 F + [/3 - a - x(2 + a + /3) ± 2 7 (1 - x 2 ))d x F 
E + 1 2 -^±P f a + P 



+ 1 ± 7(/3 - a) T ix{a + /3 + 2 ± 2s) 



F = 0. 



(46) 



The above differential equation is closely related to that of the Jacobi polynomial given by the Rodgrigues formula 

C— 1 ) n 

Pt 0) {x) = - x)~ a {l + x)-^ [(1 - xT- n {\ + xf +n ] , (47) 

which satisfies 

[(1 - x 2 )d 2 x + [/3 - a - x{a + + 2)]d x + n(n + 1 + a + 0)] p£»-« = 0. (48) 
Expanding the functions U and V as infinite series of Jacobi polynomials, 



r=0 

00 



(49) 



r=0 
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and using the recurrence relations satisfied by the Jacobi polynomials, we get the following two recurrence relations 
for A^ and B^ (where the upper sign is used if G = A and the lower if G = B): 



1 2 \ 2 ) a + (3 + 2 



G*°>± 



±- 



4 7 (r + a + l)(r + /3 + l)(r + l=FS+^ 
(2r + a + /? + 2)(2r + a + /? + 3) 



£■ + 7 2 - ( + ^rr- I 



G (r+1) -p 



47(a+l)C8 + l)(s±5 + l T . 
(a + /3 + 2)(a + /3 + 3) 



4 7 r(r + a + /3)(V±s + ^ 



(2r + a + l)(2r + a + /3) 

2 73 (a-/3)(a + ^) 
(2r + a + /?)(2r + a + /3 + 2) 



G (1) = 0, (50) 

G (r-1) 

G (r) = 0. (51) 



Determination of s Ej m 

For the series in Jacobi polynomial to converge, the constant E should satisfy a transcendental equation. To enforce 
this, we define the following quantities: 



(r) _ 2 1 {r + a){r + P){2r + a + P-2s) A^ 
(2r + a + P)(2r + a + /3 + l) A^- 1 ) 



2 7 s(a - j3)(a + , 



{2r + a + (l){2r + a + 13 + 2) 
(r) 4r 7 2 (r + a)(r + /3)(r + a + /3)(2r + a + /3 + 2s)(2r + a + - 2s) 
(2r + a + /3) 2 (2r + a + /? + l)(2r + a + /3 - 1) 

where none of the A^ r 's vanish. Eqs. (|50|) and (|5T|l may be written as 

N W _ A -(o) + ^ = o 

and 

L(r) 



jy(r+l) = K (r) _ E + 



which by iterating gives us 

N (r+1) = K (r) _ E + 

Equation ([55)1 can also be written as 
which would then give us 

N (r+l) = 



L (r) 



K (r-1) _ E . 



i('-i) 



K (r-2)_ E+ _ 



I>-2) 



k(°) -E 



L(r+2) 



I> + 3) 



Equating Eqs. (|57|) and ([5T))) . we get a transcendental equation for namely 

L ir) 



E =K {r 



K (r-l) _ E 



K' r_J J — E-\ h 



Jf(0)_E 



K(r+1) _ E 



K (r+3)_ E + . 



(52) 
(53) 
(54) 

(55) 
(56) 

(57) 



(58) 



(59) 



(60) 
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We solve this in Mathcmatica to machine precision. Convergence of the infinite continued fraction in Eq. (|59|) is 
treated in the references mentioned in [24[ ■ 

Normalization of S(x) 

We have already expanded the S(x) as a sum over Jacobi polynomials and have the recurrence relations for the 
coefficients in the expansion where all the coefficients depend on either B^ or in Eq. (j49]). Now we have, 

l-x\ a/2 fl + xV /2 ^ 

r=0 



r=0 



For the eigenfunctions to be normalized, we need two equations for the two unknowns, and B^°\ We get one 
from the Eqs. ([4"4")l and P5|) which gives us 

V = e^U. (62) 

Using Eq. (|49f at .t = 1, we have 

I ^ AW Ha! J 7 I ^ r!a! J ' 1 ' 



The other equation is obtained by using 

S 2 dx = — (64) 

Z7T 

We use (2tt)~ 1 here so that S{x)e lm ^ is normalized to 1. Using this normalization condition along with Eq. (|6Tj) . we 
have 



^ m £ ^>m I , (^) (^) ^ - 5. <«) 

which gives us 



2^ 4(0) 5(0)271 + a + ,5 + 1 n!I> + a + /3 + l) ' ~ 2n [ ' 



Hence, the second equation is 



1_ _ A jWff 2 r(n + a + l)r(n + /j + l) 

IR(0) - 27r Z^ 4(0) R(0) 



,4(0)^(0) ^(o) ^(0) 2n + a + B + 1 nllYn + a + (3 + 1) 



The ratios A^/A^ and B^/B^ are easily calculated by using their recurrence relations Eqs. (|50"|) and ([5T]) and 
using A( r > = B( r > = for r < 0. Therefore Eqs. (|63l) and (|67|) correctly determine the first coefficients in the expansion 
by choosing the A' ' whose real part is positive. 

2. Spin-weighted spheroidal harmonics as a sum over s Yi m 

The spectral decomposition of spin-weighted spheroidal harmonics in terms of spin-weighted spherical harmonics 
as in (25[, [26| has the form 

oo 

S(9)e im *= £ b js Y^ m {6,4>), (68) 
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where £ m ; n = max(|m|, |s|). Substituting the above in Eq. (|37|) . we have 

b 3 [i 2 cos 2 9 - 2 7 s cos 9 - j(j + 1)] s Y j<m = -E ^ bj S Y 3 . 



(69) 



We now write cos" 9 for n = 0, 1, 2 as sums of spherical harmonics. It is then easy to group them and the spin- weighted 
spherical harmonics as follows: From the relations 



2\/ 3^1,0, cos^ 



4V^F 1 



and 

we find (for n = 0, 1, 2) 
cos™ 6» s F^ m = 



2£+l 
An 



2n 



2^ + 1 ^' _s ' n > l-?' -s)(^,m;n,0|j,m) s Y,, m + s Y gm + 5 nfis Y im . 



(70) 
(71) 

(72) 

(73) 



j=l-n 

Multiplying Eq ([69]) with s Yk,rm integrating over the 2-sphere, and using Eq (|73|) we have 

bk-2 [l 2 Ck-2,k,2\ + b k -i [l 2 c k -iM,2 - 2sjCk-i,k,i] + b k [l 2 c k ,k,2 - 2sjc k ,k,i - k{k + 1)] 
+bk+i [l 2 c k +\.k,2 - 2.57^+1^4] + b k+ 2 \l 2 Ck+2.k.2\ = -Eb k , 



(74) 



where 



Cfc, 7, 2 



4j 
3 

Cfc, 7,1 



2j + l 



3 V 2fc + 1 



(i,m; 2, 0|fc,m)(i,-s; 2, 0|fc,-«) 



2j + l 
2fc + 1 



0', m; 1, 0|fc, m) (j, -s; 1, 0|fc, -s). 



(75) 



Eq. (|74[) can be written as a matrix equation where the b are the matrix's eigenvector and E are the eigenvalues. It 
is then easy to solve the matrix equation for the eigenvectors and eigenvalues as the matrix is band diagonal. Here 
(ji, mi; J2, m2\j, m) are the Clebsch-Gordan coefficients. 
To calculate s Yi m to high precision, we used the following analytical forms of spin-weighted harmonics at 9 = tt/2. 

1 £ + m even 



Introducing the symbol eg r , 



0, I + m odd , 



we can write 



l%m(?,0) 



2^ m (-,0) 



u+m)/2 / (2l+l7 ^(£-m)\(£ + m)\ 
K ' V 4tt U-m)\\U + m)\ 



,_ u (t +m )/2 / (2l+l)(l-m)!(l + m)! 
[ ' V 4tt£(£+1) 

f _ n «+m)/2 l (2£+l)(£-m)\(l + m)\ 
[ ' \j Att{£ -!)£{£ +!){£ + 2) 



me^m 



1 ee,m+i 



(£ - m)\\{£ + m)\\ (£-m- 1)U(£ + m- 1)!! 
[2m 2 - £(£ + 1)] ei t m 2ime£, m +i 



(£ - m)\\{£ + m)U {£ - m - l)\\{t + m - 1)!! 



(76) 

(77) 
.(78) 



Eq. (|76p is quickly obtained from the corresponding equation for P™(0), given, for example, in Arfken and Weber 
[271 ]; the corresponding relations ([77]) and ([78]) for spin- weighted harmonics follow from their definition (|40"]) . with 
recurrence relations of associated Legendre polynomials used to eliminate 9 derivatives. 

Values of the angular harmonics 2Si m u at (9, <p) = (7r/2, 0), computed as a sum of Jacobi polynomials and as a sum 
over spin-weighted spherical harmonics iYt, m are listed in Table I to show the accuracy of our calculation of these 
angular cigenfunctions. 
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t 


m 


a/M 


r /M 


o„ /p( a >/5h 

sJ£,m \-ir ) 


sSl,m (sVf?,m) 


6 


3 


0.90 


2.321 


0.25240458701173892108 


0.25240458701173889664 


10 


8 


-0.80 


8.432 


-0.077523625602031470364 


-0.077523625602031470355 


15 


-14 


0.56 


3.994 


-0.39964402300714286677 


-0.39964402300714286644 


20 


20 


0.95 


1.938 


0.53866543681715165119 


0.53866543680910100221 


25 


-24 


0.75 


3.159 


-0.38677474564628361398 


-0.3867747456462831 1527 


30 


1 


0.69 


3.439 


0.018864751317113632621 


0.018864751317113632621 


35 


-29 


0.43 


4.502 


-0.12392842343512166756 


-0.12392842343512166731 


41 


38 


0.85 


2.633 


0.41329611968515525721 


0.41329611968514376716 


45 


43 


-0.42 


7.315 


-0.44047075216769495176 


-0.44047075216769495195 


50 


-47 


0.50 


4.234 


0.37391460514301075256 


0.37391460514301011670 


54 


43 


0.29 


5.015 


0.39632284051223687540 


0.39632284051223687540 


60 


58 


0.81 


2.860 


-0.46730374640820672293 


-0.46730374640803870557 


65 


42 


-0.40 


7.255 


-0.35895613544467811490 


-0.35895613544467811490 


70 


-70 


0.67 


3.529 


-0.83907652666117403240 


-0.83907652666117381586 


75 


74 


0.80 


2.910 


0.32010913665973882714 


0.32010913665976780184 


80 


78 


0.64 


3.660 


-0.54551806835154401208 


-0.54551806835154351539 


85 


-85 


0.55 


5.555 


-0.88887357826499020752 


-0.88887357826499020734 



TABLE I: For each listed value of £, m, a and ro, we give the value of 2S^ m , obtained by using the formalism given in 

subsections C 1|) and C 2|) . with 7 = amQ = ar ^, i 2 I 1 and 6 = n/2. The fractional accuracy increases with increasing ro, 

V + a 

and all except the last five values of ro are chosen to be within a few percent of the innermost stable circular orbit for a given 
a. 



D. Numerical methods - radial harmonics 

1. Teukolsky equation 

We integrate Rh and Roq from the horizon and infinity, respectively. The homogenous solutions Rh and Roo at 
the horizon and infinity are given by the following series 

Rh - V^H^ ' (79) 

00 ? 

Roo = e iur * V 7 ; " tt , (80) 

where 

r\+ a 2 r 2 _ + a 2 
r+ = r H — — In r — r+ In r — r_ and (81) 



r± = M ± VM 2 - a 2 . (82) 
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The expansion coefficients satisfy the following recurrence relations 



c n =< 2i(n — 6)uMc„_3 + 



E - n 2 + 7n - 12 + uj 2 a 2 + i{8n - 3A)uM^l - a 2 /M 2 + i(An - 1A)wM 



C n -2 



Amua + (2E - An 2 + 18n - 18 + 2wV) A /l - a 2 /M 



+ H2(n - 2)uiM - iA(2n - 5)uia 2 /M - iAma/M + i\2(n - 2)uMy/l - a 2 /M 2 



Cn-l 



(An 2 - 8n)(l - a 2 /M 2 ) + m 2 a 2 /M 2 - Amwa - Amu ay/ 1 - a 2 /M 2 



i[8nMu)(l - a 2 /M 2 ) + {8nujM - Ama/M)y/l - a 2 /M 2 



i(n+l)a 4 i(n + l)(2n + l + iuoa 2 /M)a 2 

2a,A/5 n ~ 5 + ^Jm~* "~ 4 

-AmMa + (An + 8)uMa 2 - i(An 2 + 8n)M 2 + i(E - 2n 2 -An- m 2 )a 2 + icv 2 a 4 

2nwM 2 

2ma - (2n + l)wa 2 + i(-E + 2n 2 + 5n + 3 + 2mwa - uj 2 a 2 )M 



d n - 



TILUM 

(An - 4)wM + i(E - n 2 - in - 2 + cu 2 a 2 ) 
2nuj 



d n -i- 



(83) 



(84) 



We use a 7th order Runge-Kutta routine to solve for the homogenous radial solutions using the above initial/boundary 
conditions, obtaining values of ?/>5 ct to an accuracy of 1 part in 10 13 . 



2. Sasaki-Nakamura equation 



The Sasaki-Nakamura equation is 



d 2 X ^dX 



dr* 2 dr 



F- UX = 



where X is related to the radial part, R4 of p "04 by 



i?4 — — 
77 



where 



AX 



X 



y/r 2 + a 2 

The radial parts, R4 and i?o are related to each other by the relation 



Ro - c^, 



where c is a constant. The function F is 



F 



where 



j) = - \2iuM + A(A + 2) - 12au(auj - m) + 



i A 

n r 2 + a 2 

8ia[3au; — X(aco — m)\ 



(85) 



(86) 



(87) 



(89) 



+ 



-24iaM(ao;-m) + 12a 2 [l - 2(aw - m) 2 ] 2Aia 3 (auj - m) - 2AMa 2 12a 4 



(90) 
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and the prime denotes a derivate with respect to r. 
The function U is 



(r 2 + a 2 ) 2 



where U\ and G are 



A 2 



— \2a + — ) «+-r 

dr \ A J r] \ A 



and 



G = - 



2(r — M) 



-A 



(r 2 +a 2 ) 2 ' 



+ a 

The a, ft and V that appear above are given by 

' K 2 + AiK(r - M) 



V=- 



ft = 2A 



A 



-iK + r - M 



+ 8iuir + A , 



2A 



-iKft . . x 6A 
A r 



where A' = (r + a )u> — am. We use the following boundary conditions at the horizon and at infinity, 



where 



and m 



n=0 



2Mr 4 



r = r 



log 



2M 



27\/r_ 



log 



2M 



(91) 



(92) 



(93) 

(94) 
(95) 

(96) 



d 
<h 

d 2 

d 3 



-i(2 + A + 2omw) 
2^ ' 

A 2 + A(4 + amu) - 12iMui + 4amu(l + amuj + 2iMuj) 

8J 2 ' 
i\ 3 i\ 2 (3amuj-l) i\(-2 - 2amuj - ZiMoj + 2a 2 ui 2 + 3a 2 m 2 uj 2 + QiamMuj 2 ) 



48w 3 24w 3 12w 3 

i[6iM - 6am - 3acj(a + imM) + 2amw 2 {a 2 - AM 2 ) + (amuj) 2 (am + 6iM)} 

6ui 3 

A 4 A 3 (8amo;-12) A 2 (12 - 72amw + ASiMu + 32a 2 cj 2 + 24a 2 m 2 uj 2 + A8iamMcu 2 ) 



384w 4 



384w 4 



384w 4 



A[80(l - amw) + 2881Mw + 128a 2 w 2 (amw - 1) + 16 (amu) 2 {2amuj + UiMui - 7) - 256amA/ 2 w 3 ] 

384w 4 

H ^[30am - 30iM - 6a 2 cj - 15cj(am) 2 + WiamMuj + 45M 2 u - l6a 3 moj 2 - 2(am) 2 uj 2 

24w 3 

- 18iM(aw) 2 - 6iM{amu) 2 + Aam(Muj) 2 + 8a 4 m 2 uj 3 + (amfuj 3 + 24imM(au) 3 + \2iM(amu) 3 

- Uoj(amMuj) 2 - A8iam{Muj) 3 ] . (97) 

The accuracy of our radial eigenfunctions is shown in Table II, which exhibits values of \Ri n R ou t/W\ computed by 
independently integrating the Tcukolsky equation and the Sasaki-Nakamua equation. 
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I 


m 


a/M 


r /M 


Teukolsky 


Sasaki-Nakamura 


56 


53 


0.50 


4.30 


3.32878980028437xl0" 2 


3.32878980028426 xlO" 2 


85 


84 


0.50 


4.25 


2.21479032560870xl0" 2 


2.21479032560832xl0" 2 


83 


76 


0.30 


5.00 


2.66435949501300xl0~ 2 


2.664359495012799 x 10 -2 


42 


1 


0.30 


6.00 


5.774432559632249 x 10~ 2 


5. 774432559632295 xlO" 2 


6 


-1 


-0.25 


6.80 


4.40438489554819xl0 _1 


4.40438489554828xl0 _1 


56 


37 


0.61 


3.79 


2.5750354691478 xlO" 2 


2.5750354691415 xl0~ 2 


65 


-58 


-0.76 


8.33 


6.068437461161814xl0~ 2 


6.068437461161780 xlO^ 2 


83 


63 


-0.71 


8.18 


4.5482452073276329 x 10~ 2 


4.5482452073276363 x 10 -2 


76 


74 


0.92 


2.20 


1.056332700247276 xl0~ 2 


1.056332700247222 xlO -2 


68 


62 


0.16 


5.50 


3.62341811077144xl0" 2 


3.62341811077135xl0" 2 


79 


67 


0.70 


3.40 


1.69778531677806xl0~ 2 


1.69778531677799xl0" 2 


56 


-42 


-0.25 


6.80 


5.4130389893688027 x 10~ 2 


5.4130389893688006xl0~ 2 


79 


71 


0.58 


3.95 


2.09651683213181xl0" 2 


2.09651683213165xl0" 2 


82 


75 


-0.2 


6.70 


3.77916963861862xl0" 2 


3.77916963861852xl0" 2 


54 


35 


0.14 


5.55 


4.31642221118258xl0" 2 


4.31642221118267xl0" 2 


67 


52 


-0.64 


8.00 


5.497676877828166 xlO -2 


5.497676877828151 xlO -2 


75 


61 


0.78 


3.05 


1.52296461688768494x 10 -2 


1.52296461688768459 xlO -2 


77 


68 


-0.14 


6.50 


3.85508542936437xl0" 2 


3.85508542936429 xlO" 2 


65 


53 


-0.72 


8.30 


5.946127670533624 xlO -2 


5.946127670533644 x 10" 2 


81 


66 


-0.58 


7.80 


4.456906403886178 x 10~ 2 


4.456906403886146 x 10~ 2 


86 


81 


0.47 


4.35 


2.20499986904938xl0" 2 


2.20499986904907xl0" 2 


53 


-53 


-0.23 


6.80 


6.057956364463913 xlO -2 


6.057956364463899 x 10" 2 


70 


65 


0.39 


4.70 


2.94548366362787xl0~ 2 


2.94548366362775 xlO" 2 


52 


33 


0.04 


6.15 


5.053852006520149 x 10" 2 


5.053852006520189 xlO -2 


44 


44 


-0.34 


7.10 


7.645731150733178 xlO -2 


7.645731150733167xl0 -2 


70 


68 


0.21 


15.00 


1.027161976019888X10" 1 


1.027161976019891 xlO -1 


72 


67 


-0.66 


43.50 


2.9619125920183531 x 10" 1 


2.9619125920183542X10" 1 


40 


-39 


0.21 


25.00 


3.02127126128480971 xlO^ 1 


3.02127126128480916xl0 _1 


61 


58 


0.91 


72.50 


5.849950541617765 xlO" 1 


5.849950541617792xl0 _1 


50 


2 


-0.88 


15.50 


1.4349668057921668X10" 1 


1.4349668057921630X10" 1 


65 


12 


-0.44 


100.00 


7.558239183584133 xlO" 1 


7.558239183584106 xlO" 1 



TABLE II: For each listed value of £, m, a and ro, we give the value of |i?i n -Rout/ W\, obtained by integrating the Teukolsky 
equation and the Sasaki-Nakamura equation. As in Table Hi C 21 the fractional accuracy increases with increasing ro, and all 
except the last six values of ro are chosen to be within a few percent of the innermost stable circular orbit for a given a. 

III. COMPUTATION OF THE PERTURBED METRIC 

We can now use Eq. (f29|) to compute the perturbed metric in an ORG in terms of the Hertz potential 4'. We first 
convert the spin-weighted spheroidal harmonics in Eq. (|28| to spin-weighted spherical harmonics as follows: Restoring 
the suppressed indices s, £, and m to S and bj in Eq. (pJSf . 



S — sSlmi bj — bgj m , 



(98) 
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with s — 2, we have 



Chi 



; E 

Cm 



L> 2 + U4M 2 m 2 n 2 
{-l) m DR tm + 12imMClRe m ^ innt 



D 2 + U4M 2 m 2 n 2 



EE 1 



{-l) m DR lm + 12imMnR em _ e _ lmnt 



D 2 + lUM 2 m 2 n 2 



(99) 



To write the nonzero tetrad components of the metric perturbation of Eq. (|29j) , we replace spin-coefficients by their 
values for Kerr given in Eq. (j3J), we write the derivatives operators A and D in the explicit forms implied by ([3]), and 
we use Eq. (|41[) to replace the angular derivative operator 5 by its form in terms of 3. We thereby obtain 

/in = ^2 [§ 2 - 2ia(g + imQ.) sin 9 3 - 3a 2 g 2 sin 2 6 + a 2 {g - imQ)(3g + imfi) sin 2 9]^ lm2 Y lm + c.c, (100) 



'13 



Cm 



E 

Cm 



-1 A 

4V2Xg 3 



2<9 r 8 — ia — g + 2im£l — 2i— cos 6* J sin#d r — 2 I g + g + 4 — — \- iq\<5 



A 

r-M 



Qiag 2 sm9 — ia (2g + g — imfl + 2i— cos O^j y3g + 4 — ^— + iqj sin 9 

( _ r-M \ 
— ia (3g + imO) g — 2g — 4 — ig sin 9 



t< m2 F <m , (101) 



= E 7 



1 A 5 



/:'m 



4E 2 g 4 



« + [q - Q 



r - M 2r 



— + 2iq )d r + 3g 2 + — - 



r-M) 2 q(r-M) + m£lr 



2g+g-A- 



r-M 



A 



2- + 2iq) [3g + 4 



2i 



r-M 



A 



^ImiYtm. (102) 



Because \& has spin- weight 2 and has spin- weight 1, 3 acting on VP and on Qvp has the form (j4Tj) with s = 2 and 
s = 1, respectively. 

We next expand the operators acting on ^ in powers of the small parameter cos# and show below that terms 
involving cos 9 and cos 2 9 do not contribute to h 1 ™ (or to the nonzero part of the singular field). Using Eq. (|40| for 
the action of 3 on 2^m, we obtain 



EM E 



'•Attn 1 



.Yl, 



(103) 



Cm 



where the radial functions s ^£ m & re given by 



/ (I + 2)! r 2 



!^° m = -V(^-l)(^ + 2)ar(i + mfir), 
2 A° hn = ^a 2 mQr(2i + mfir), 



„i • /(* + 2)l 

o-Afm = -HI— — — ar > 



(1-2) 



iA\ m = -^(l-l)(e + 2)a 2 {l - 2imQr), 
2 A\ m = mf2a 3 (l — imrQ), 



2 _1 (£ + 2)! 
fm 2V(^-2)! 



= 1^-1)^ + 2)0^ + ™U(r 2 + 2a 2 )], 
2 A 2 m = -lmfia 2 [2zr + mfi(r 2 + a 2 )]; 



2 2 
£m s— 1 n— 

where, with <? defined by Eq. (|14p, the functions S S|L are given by 



A = -J ie l) ^ + 2) [rA^ m + (2r 2 - 2a 2 - iqrA)* tm ) , 



2 B° (m = \am^[rA% m + (2r 2 - 2a 2 - iqrA)* em ], 



lBL V l) ^ + 2) - r [rA% m + (2r 2 - 2a 2 - iqrA)^ m ], 

2 B\ m = (2 + 3imflr)A% m + [8(r - M) + 3mfigrA - 2ima + AimVl{2r 2 - a 2 )] V lm 



A = ^^^ ^J t 2 ^™ + (3r 2 + 2Mr - 5a 2 - 2iqrA)* lm ] , 

2 Bj m = _LiL{[ TO f7 r ( r 2 + 8a 2^ + i2ia 2 ]Aty' em + [12ma 3 + 2mO(r 4 - a 2 r 2 + 4Ma 2 r - 16a 4 )]* £m 

+ t[48a 2 (r - M) - mClqr(r 2 + 8a 2 )A]* fm }; 



and 



h 33 = J2J2 C ^n^s n 9 2 Y tm , 



with C2 m given by 



C\ m = ±A 2 ^ m + ~(3r 2 -2Mr-a 2 -iqrA)V hn 



r 2 - 2M 2 -a 2 + 2— - -q 2 A 2 - *^(2r 2 - Mr - a 2 ) 
r 4 2 r 



18 



r 1 



-iaA 2 ,, iaA . 

v'L * 



ia 

$i 

r 



3A 

M) — 2im\a - Q(r 2 + a 2 )} 

2r 



2imnAr + 8(M - 2) 2 - 2A + + 6im[a - n(r 2 + a 2 )](M - r) 

r 



3Aim[a - n(r 2 + a 2 )} 



[o-n( 



2 i 2 

r + a 



(109) 



-2a 2 A 2 „ Aa 2 , 
.iff'! vp' 



r 



K A 

16(r - M) 4im[a - fi(r 2 + a 2 )] 

2r 



4imfiAr + 16(M - r) 2 - 2A + 10MA + \2im{M -r)[a- Q(r 2 + a 2 )] 



5A[a - Q(r 2 + a 2 )} 
2~r 



2m 2 [a- Q{r 2 + a 2 )} 2 



(110) 



We now argue that the terms involving cos™ 9 with n ^ in Eqs. (|103[) . (|105[) and p07p can be ignored. We refer 
here to our description in Sec. Ill Bl of the the renormalization procedure in terms of ip® and VP. The components 
/i™ 1 of the renormalized radiation-gauge metric are given by Eqs. ()103[) . ()105[) and (|10T[) . with ^ replaced by \l/ ren . 
Because vl/ rcn is smooth and hence finite at the particle, no term proportional to cos# contributes to h 1 *™. 

It can happen, however, that the cos# terms contribute to h^ v , written in terms of The cos 2 9 terms cannot 
contribute, because they are 0(p 2 ) and multiply terms whose sum is at most 0(p^ 1 ) (terms involving two derivatives 
of fy). Similarly, the cos# terms cannot contribute to the leading term in h , because they are one order in p smaller 
than the leading term in H s . At subleading order, however, they give an O(p ) contribution that has odd parity. 
The parity of this contribution to the singular field follows from that fact that, at leading order, \I' ret is even under 
parity about the position of the particle (shown in detail in Sec. HID of [l|), while cos 9 is odd. The leading-order 
contributions to the cos# terms come from terms involving two derivatives of <J> rot , and these are again even under 
parity, implying that at leading order, the contribution from r > cancels the contribution with opposite sign from 
r < r . Finally, lower-order contributions multiplying cos# arc order p, vanishing at the particle. 

One can then compute H ren by subtracting the leading part of the singular field (which coincides the Lorenz-gauge 
singular field) and by omitting terms that involve cos 9. We have verified this agreement numerically (see Sec. IVII) . 
finding that the order L° part of Hf agrees to one part in 10 12 with its analytic form computed by Linz [23| and that 
the O^L" 1 ) contribution vanishes to within the accuracy of the computation. 



IV. GAUGE-INVARIANT QUANTITIES 

In this section, we obtain expressions for the related quantities AU and Afi that give, respectively, the change 
in the redshift factor of a trajectory at fixed angular velocity and the change in the angular velocity of a trajectory 
at fixed redshift factor U . Each of these quantities is invariant under gauge-transformations generated by helically 
symmetric gauge vectors £ Q and each can be written in terms of the similarly gauge-invariant quantity 

H ren := ^h%$u a u p , (111) 

where h™p is the renormalized metric perturbation. 

As shown by Mino et al. [l9| (see also Quinn and Wald [2(| and Detweiler & Whiting [IH), at order m/M the 
particle moves along a geodesic of the metric g a p + h r ®g, where g a p is the background (Kerr) metric. Denote by 
u a = U(t a + il<f) a ) the particle's 4- velocity, normalized with respect to g a p + h™p, 

{g a p+h^)u a u p = 1. (112) 

We consider first the difference AU between the value of U for a circular geodesic of the perturbed metric and its 
value at the circular geodesic of the unperturbed metric with the same value of angular velocity f2. Formally, because 
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perturbations to first-order in ra/M are linear in m, we can write 

AC/ = mJU/(m,Q)| m = . (113) 

We will denote by 5U the gauge-dependent change in U at a fixed value of r. Denoting by U(m,r) value of U for 
the circular geodesic at radial coordinate r of the metric g a p + h™, we have 

d 

5U:=m— C/(m,r)| ro=0 . (114) 

Let £ Q be the radial vector joining the unperturbed circular geodesic with angular velocity fl to the perturbed geodesic 
with the same angular velocity: Formally 

fi(m = 0, r) = fi(m, r + O + 0(m 2 ). (115) 

Then 

AU = SU + £$U, Afl = Sn + £ € Q = 0, (116) 

where AJ7 and <50 are defined as in Eqs. (|113[) and (| 1 14|) . 

We can now quickly compute AU from Eq. (|112p . showing as follows the relation 

AU = -u*ir cn . (117) 

Define k a in the equatorial plane by k a = t a + £l(f> a , with fi = il(m, r), and let k a be the Killing vector k a = t a + ilo4> a , 
where O is the angular velocity of the unperturbed orbit through r = r . Then, applying A to the normalization 
equation (|112[) and evaluating the expression at r — r , we have 

= A(g a0 U 2 k a kP) = {S + L^ig^U 2 ^) = {h™$ + £^g af i)u a u^ + ^AU + 2Uu a (S + £ € )fc Q , (118) 

with S again the perturbation at fixed radius r, as in Eq. (|114[) . Using the fact that is helically symmetric, we 
now see as follows that the terms £^g a /su a u 13 and (5 + £^)k a vanish. At r = r , we have U(m = 0,r ) = it*, whence 
u a = u f k a , and Eq. ([33)) then implies £^g a /su a u 13 = 0. Finally, because the coordinate system is independent of m, 
we have 5t a = = 8(j> a (that is, dt and do not change), and the last term vanishes: 

(5 + £ € )fc Q | r =r = £it a + n £^ a = £$k a = -£ k C = 0. (119) 



1 



,t h ren a 



'a/3 



U W 



From the two surviving terms on the right of Eq. (|118[) , we obtain the claimed form AU = 

The change in the angular velocity at fixed U is similarly gauge invariant and is easily obtained from AU. With Q 
regarded as a function of m and U we define its change at fixed U by 

Att:=m^(m,[/)| m=0 . (120) 

From the fact that, at fixed m, S7(m, U) is the inverse of C/(m, f2), it follows that 

d 

AO = — — :fi(m = 0, U)AU, (121) 
a U 



Aft = l — H Icn . (122) 



implying 



The resulting values of AU and AO are presented in Tables [TTT1 and HVl 

For completeness, we give here explicit expressions for the quantities SU and SD,. These depend on the gauge- 
dependent acceleration a a , the self- force per unit mass, which is ordinarily defined with the perturbed trajectory 
parametrized by proper time with respect to the background metric, implying for the 4- velocity u a the normalization 



g a p u a uf = l. (123) 
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Denoting by r and f proper time along a trajectory with respect to g a p and g a p + h™g, respectively, we have 
df 

u a = — u a , with 
dr 

— = 1 - -h^u a uP = 1 - H Icn . (124) 
df 2 p 

The acceleration of the perturbed trajectory with respect to the background metric g a p is given by 

o 

a a ■= u^VfjU 01 (125a) 
= -(g aS + u a u 5 ){V 1 h P 5-\v s hp 1 )u f) u 1 , (125b) 

o 

where V a is the covariant derivative operator of g a p. 

Using Eqs. (|123[) and (|125a[) . we find for the changes in f2 and u* of a trajectory at fixed radius r 

8n = n 2aMl ' 2r \+f^- iM \ r (126) 

2M(rl /2 + oAfVa) 

and 

Su l = u l V \ - °-^a r , 127 

where Q and u are given by the relations in Eq (j6]), f2 is the frequency measured by the observer at infinity. Here 
one is comparing the values of and it* for circular geodesies of the perturbed metric g a g + to their values for 
a circular geodesic of the unperturbed metric g a p at the same value r = tq of the radial coordinate. The expressions 
are valid in any gauge, but the values of a r , 60, and 6U are gauge-dependent. 



V. LOWER MULTIPOLES 



The metric recovered from ipQ en specifies the perturbation up to the contribution that comes from the change in 
mass and angular momentum of a Kerr metric and from a change in the center of mass that is pure gauge except at 
r = r (and that does not contribute to H lcn ). 2 In this section we calculate the contribution to the gauge-invariant 
fjret f rom change in mass and angular momentum due to the presence of the orbiting particle of mass m. We 
calculate them in the "Kerr gauge" ; that is, they are written as the first-order perturbations of the Kerr metric in 
Boyer-Lindquist coordinates associated with the changes 5M and 6 J in its mass and angular momentum. These two 
parts of the metric perturbation are thus stationary and axisymmetric, and they are associated with a stationary, 
axisymmetric part of the stress-energy tensor. For a particle in circular orbit, SM and 6 J have the simple forms 

SM = E = mu a t a , 6 J = L = -mu a (j) a (128) 

(for our H signature), as stated by L. Price (28|. The expression for S J follows, for example, from the Ko- 

mar formula for angular momentum, valid for a stationary axisymmetric perturbation of a stationary axisymmetric 
spacetime; it implies 

6J = - [ T^^dSa = -mu^, (129) 
Jv 

when there is no change in the angular momentum of the black hole. The change in the mass follows from the 
Bardeen-Carter-Hawking first law of thermodynamics for black holes and matter, which gives 

SM = nSJ + m/u* = mu t . (130) 



2 In the mode-sum renormalization, the individual modes of the metric are computed as the limits of their values as r — > ro from r < tq 
or r > ro. Because the because the metric perturbation associated with a change in the center of mass is pure gauge for r ^ ro, these 
limits vanish. 
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(One can also use the generalization of the first law to helically symmetric binaries [29|-|3lJ , but the stationarity and 
axisymmetry of the relevant perturbation again means that the original form of the first law is valid.) 

Although one can also find SM directly from the Komar expressions in terms of the timelikc Killing vector shared 
by the background and perturbed spacetimes, a warning is needed: An unexpected subtlety arises in using the Komar 
expression for the mass. The Komar mass on a sphere outside the particle orbit has the form 

SM = 5-^-1 V a t' 3 dS aP , (131) 

47T J 

where dS a p = ^eap^sdS"^ (implying dSt r = ^\/—gd9d(f>). If one requires that the change in the Komar mass of the 
black hole vanish and assumes that St a = 0, the result is an incorrect expression for SM. Gauss's theorem gives 

SM = 5^-1 V f3 V a t p dS a + 5^- I V a t^dS aP 

47r JV 4?r ./horizon 

= f {2T a p ~Sf 3 T)t f} dS a + 5^ f V a t^dS aP , (132) 

J J horizon 

and 

J (2T a p - S%T)tPdS a = m{2u t - l/u*) ^ mu t . (133) 

The discrepancy arises from a rescaling of h u = h a pt a t^ near the horizon for a time-independent perturbation. This 
is easiest to see for a Schwarzschild background, where the perturbation in the mass arises from the spherically 
symmetric part of the perturbation - from the perturbation due to a spherical shell of dust whose particles have 
trajectories isotropically distributed over all circular geodesies at r = r$. The Komar mass is gauge invariant, 3 and 
we can can compute it in a Schwarzschild gauge. The perturbed field equation then requires continuity of htt ='■ e 2 * 
across r = r , and 5$ is constant inside r — r , implying a constant rescaling of time for r < r . The result is that 
the expression for the change in the Komar mass at the horizon is evaluated with a rescaled metric but without a 
rescaled t a , giving a nonzero result, 

5^- [ V a t?dS al3 = m(l/u' ~ ut), (134) 

-/horizon 

that yields the correct value SM = mut for the change in the spacetime mass. The change in h a /3 inside r = has 
the form £%g a p for a vector £ Q linear in t; because t a remains fixed, however, this is not a gauge transformation of 
the integrand y a t^dS a p (the integrand does not change by £^(V Q t /5 'dS a p)). 

Finally, Eq. (|129j) for the change in angular momentum is valid, because the rescaling does not alter the Komar 
expression for the angular momentum at the horizon. 

To calculate the metric perturbation that comes from the change in mass and angular momentum, we find the first 
order perturbation of the (relevant components of the) Kerr metric in Boyer-Lindquist coordinates which are 



h 



2SM 



ii 



ht<t> = 

2(M + r)aHM 

^ = M~r ( 35) 



for the change in mass, and 



ht<p 





2MSJ 
r 

2a{2M + r)SJ 
Mr 



(136) 



3 Gauge invariance of the integral (11311) . over a sphere S where there is no matter, can be seen as follows. Let S' be another sphere 
homologous to S with no matter in the region between them. Then, for any given choice of gauge, the value of the Komar integral is 
the same on the two spheres. Consider a gauge transformation associated with an arbitrary gauge vector £ Q defined in a neighborhood 
of S, and extend smoothly so that it vanishes on 5'. In the new gauge, the value of the Komar integral on S' has not changed; and 
it must again have the same value on S and S". Its value on S is therefore unchanged. 
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for the change in angular momentum. From these expressions and Eqs. (|6|) for u* and f2 (and using = fiu 4 ), we 
obtain 

_ m(r 2 ± 2oMVV^ - a 2 )(r 3 / 2 - 2Mr 1 / 2 ± aM 1 / 2 ) 

HSM ~ r 9/4( r 3/2 _ 3Mr l/2 ± 2aM 1 /2)3/2 ( 137 ) 



and 



M i/2 m(r 2 -p 2qM 1 /V/ 2 + a 2 )(±q - 2M 1 / 2 r 1 / 2 ) 

" - r 9/4( r 3/2 _ 3 JV/rV2 ± 2aM 1 /2)3/2 ( ) 

where the upper (lower) sign is used for direct (retrograde) orbits. 

VI. NUMERICAL RESULTS 

The renormalization of H follows the mode-sum method described in Sec. Ill Bl After the odd-parity terms - terms 
involving cos# in Eqs. (|103[) . (|105[) and (|107[) - are omitted, the method is identical to that used in [J] for a particle 
in circular orbit in a Schwarzschild background. 

In Eq. with the odd-parity part of H gone, the remaining 0(L~ 2 ) terms vanish at the particle, allowing us 
to write the remaining part Hf of Hf in the form 

'■^gw (139) 

where p2k(&) is a polynomial in I of order 2k for which 



OG 1 



We numerically match Hg (where the tilde again denotes a value computed with odd-parity terms omitted) to this 
expansion of Hf, 

k 

H s = Eo + y^k (141) 

and extract the regularization coefficients Ei^ up to fc max between 8 and 10. The method used in the numerical 
matching and an error- minimization criterion for the choice of fc max are described in detail in Q ■ The resulting value 
of H len is given by 



^rcn = g [jff rct _ = g [tfret _ £|] ) (142) 



£=0 fcO 

74. The analytical value of Eq [23| is given by 



-^analytical — F =-K" f : 7 ) , (143) 

Try/(1 + P)g gg \1 + PJ 

where (3 = 9 ** ^" +L and K is the complete elliptic integral of the first kind, K(m) = (1 — msin 2 4>)~ 1 / 2 d4>. 

We compare the value of Eq obtained by numerical matching to the above analytical result and observe that they 
agree to 12 significant figures. 
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ro/M 


a = -0.9M 


a = -0.7M 


a = -0.5M 


a = 0.0M 


a = 0.5M 


a = 0.7Af 


a = 0.9M 


4 


- 


- 


- 


- 


- 


-0.39639405 


-0.32811192 


5 


- 


- 


- 


- 


-0.31443977 


-0.27861234 


-0.25156061 


6 


- 


- 


- 


-0.29602751 


-0.23463184 


-0.21756347 


-0.20361838 


7 


- 


- 


- 


-0.22084753 


-0.18875155 


-0.17902998 


-0.17073001 


8 


- 


- 


-0.20415909 


-0.17771974 


-0.15838853 


-0.15222199 


-0.14680897 


10 


-0.15129436 


-0.14557511 


-0.14033900 


-0.12912227 


-0.12019572 


-0.11717475 


-0.11443451 


15 


-0.083291764 


-0.081933637 


-0.080646922 


-0.077725319 


-0.075195106 


-0.074284771 


-0.073429473 


20 


-0.058142984 


-0.057590366 


-0.057059948 


-0.055827719 


-0.054723506 


-0.054316065 


-0.053927537 


30 


-0.036504919 


-0.036334869 


-0.036169772 


-0.035778314 


-0.035416550 


-0.035279964 


-0.035147937 


50 


-0.021026283 


-0.020984416 


-0.020943414 


-0.020844656 


-0.020751199 


-0.020715285 


-0.02068020 


70 


-0.014784459 


-0.014767331 


-0.014750491 


-0.014709646 


-0.014670583 


-0.014655454 


-0.014640606 


100 


-0.010234918 


-0.010228170 


-0.010221515 


-0.010205282 


-0.010189625 


-0.010183523 


-0.010177512 



TABLE III: This table presents the numerical values of AU for different values of ro /M and a. They are accurate to a fractional 
difference of order 10~ 8 . 



ro/M 


a = -0.9M 


a = -0.7M 


a = -0.5M 


a = 0.0M 


a = 0.5M 


a = 0.7M 


a = 0.9M 


4 












0.054267340 


0.052559297 


5 










0.047924050 


0.046963137 


0.046465534 


6 








0.042727891 


0.040850942 


0.040470951 


0.040275775 


7 








0.036056740 


0.035175043 


0.034994056 


0.034900960 


8 






-0.031876878 


0.031046361 


0.030576073 


0.030478282 


0.030427413 


10 


-0.024543706 


-0.024365158 


-0.024209291 


0.023913779 


0.023742658 


0.023706031 


0.023686199 


15 


-0.014462899 


-0.014434162 


-0.014408846 


0.014359915 


0.014330238 


0.014323381 


0.014319216 


20 


-0.0098123143 


-0.0098041785 


-0.0097969635 


0.0097828022 


0.0097738694 


0.0097716692 


0.0097702125 


30 


-0.0055825069 


-0.0055810957 


-0.0055798309 


0.0055772872 


0.0055755838 


0.0055751254 


0.0055747889 


50 


-0.0026871991 


-0.0026870390 


-0.0026868933 


0.0026865907 


0.0026863727 


0.0026863081 


0.0026862560 


70 


-0.0016464854 


-0.0016464466 


-0.0016464110 


0.0016463355 


0.0016462787 


0.0016462609 


0.0016462460 


100 


-0.00097498493 


-0.00097497623 


-0.00097496816 


0.00097495060 


0.00097493692 


0.00097493244 


0.00097492852 



TABLE IV: Numerical values of MAO, for different values of ro/M and a. The values are accurate to a fractional difference of 
order 10~ 8 . 



VII. DISCUSSION AND FUTURE WORK 



The results here are based on the computation of the invariant H lcn 7 and work now underway with A. Le Tiec shows 
that AU from the computations in our modified radiation gauge agrees with the post-Newtonian series for AU linear 
in the spin parameter a/A/: A preliminary matching shows that the first coefficient in the pN series agrees to five 
significant digits. The results here also agree with those of a separate EMRI computation by Dolan, who works in a 
Lorenz gauge using an effective source method (agreement is within their numerical error bars of order 10 -2 ). Finally, 
we have also begun work to extend the computation reported here to find the self-force on a particle in circular orbit 
in a Kerr background. 
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